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On a universal photonic tunnelling time 
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We consider photonic tunnelling through evanescent regions and obtain general analytic expressions 
for the transit (phase) time r (in the opaque barrier limit) in order to study the recently proposed 
"universality" property according to which r is given by the reciprocal of the photon frequency. 
We consider different physical phenomena (corresponding to performed experiments) and show that 
such a property is only an approximation. In particular we find that the "correction" factor is a 
constant term for total internal reflection and quarter-wave photonic bandgap, while it is frequency- 
dependent in the case of undersized waveguide and distributed Bragg reflector. The comparison of 
our predictions with the experimental results shows quite a good agreement with observations and 
reveals the range of applicability of the approximated "universality" property. 



I. INTRODUCTION 

In recent times, some photonic experiments deal- 
ing with evanescent mode propagation have drawn some 
attention because of their intriguing results. All such ex- 
periments have measured the time required for the light 
to travel through a region in which only evanescent prop- 
agation occurs, according to classical Maxwell electrody- 
namic. If certain conditions are fulfilled (i.e. in the limit 
of opaque barriers), the obtained transit times are usu- 
ally shorter than the corresponding ones for real (not 
evanescent) propagation through the same region. Due 
to the experimental setups, this has been correctly in- 
terpreted in terms of group velocities Q greater than c 
inside the considered region. Although there has been 
some confusion in the scientific community, leading also 
to several different definitions of the transit time Q , these 
results are not at odds with Einstein causality since, ac- 
cording to Sommerfeld and Brillouin the front ve- 
locity rather than the group velocity is relevant for this. 
Waves which are solutions of the Maxwell equations al- 
ways travel in vacuum with a front velocity equal to c 
while, in certain conditions, their phase and group veloc- 
ities can be different from c |]lO|. It is worthwhile to ob- 
serve that the quoted experiments are carried out study- 
ing different phenomena (undersized waveguide, photonic 
bandgap, total internal reflection) and exploring different 
frequency ranges (from optical to microwave region) . 
The interest in such experiments is driven by the fact 
that evanescent mode propagation through a given region 
can be viewed as a photonic tunnelling effect through a 
"potential" barrier in that region. This has been shown, 
for example, in Ref. |^| using the formal analogy be- 
tween the (classical) He mholtz wave equation and the 



(quantum mechanical) Schrodinger equation (see also 
Ref. [Q). In this respect, the photonic experiments are 
very useful to study the question of tunnelling times, 
since experiments involving charged particle (e.g. elec- 
trons) are not yet sensible enough to measure transit 
times due to some technical difffculties [p^ . 
From an experimental point of view, the transit time t 
for a wave-packet propagating through a given region is 
measured as the interval between the arrival times of the 
signal envelope at the two ends of that region whose dis- 
tance is D. In general, if the wave-packet has a group 
velocity Vg, this means that r = D/vg. Since Vg — duj/dk 
(k wave- vector, to angular frequency), then we can write 



duj 



(1) 



where d(j) — D dk is the phase difference acquired by the 
packet in the considered region. The above argument 
works as well for matter particles in quantum mechanics, 
changing the role of angular frequency and wave-vector 
into the corresponding ones of energy and momentum 
through the Planck - de Broglie relations. 
However, difficulties arise when we deal with tunnelling 
times, since inside a barrier region the wave-vector (or 
the momentum) is imaginary, and hence no group veloc- 
ity can be defined. As a matter of fact, different defini- 
tions of tunnelling time exist. While we refer the read to 
the quoted literature [^, here we use the simple defini- 
tion of phase time which coincides with Eq. (^. In fact, 
although Vg seems meaningless in this case, nevertheless 
Eq. (|l|) is meaningful also for evanescent propagation. 
The adopted point of view takes advantage of the fact 
that experimental results |l|- seem to confirm the def- 
inition of phase time for the tunnelling transit time. 
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FIG. 1. A barrier potential V{z) for a particle or a barrier 
refractive index n{z) for an electromagnetic wave. 

Recently, Haibel and Nimtz |^ have noted that, regard- 
less of the different phenomena studied, all experiments 
have measured photonic tunnelling times which are ap- 
proximately equal to the reciprocal of the frequency of 
the radiation used in the given experiment. Such a "uni- 
versal" behaviour is quite remarkable in view of the fact 
that, although photonic barrier traversal takes place in 
all the quoted experiments, nevertheless the boundary 
conditions are peculiar of each experiment. 
In the present paper we carefully study the proposed uni- 
versality starting from a common feature of tunnelling 
phenomena and, in the following section, derive a gen- 
eral expression for the transit (phase) time. Different ex- 
periments manifest themselves into different dispersion 
relations for the barrier region. We then analyze each 
peculiar experiment in Sects. III,IV,V and compare the- 
oretical predictions with experimental observations. Fi- 
nally, in Sect. VI, we discuss our results and give conclu- 
sions. 

Note that, differently from other possible analysis (see, 
for example, the comparison with a photonic bandgap 
experiment in |p5[), we deal with only tunnelling times, 
which have been directly observed, and not with veloci- 
ties which, in the present case, are derived from transit 
times. 



II. PHASE TIME AND DISPERSION RELATION 



In this paper we study one-dimensional problems or, 
more in general, phenomena in which evanescent prop- 
agation takes place along one direction, say z. Let us 
then consider a particle or a wave-packet moving along 
the z-axis entering in a region [0, a] with a potential bar- 
rier V{z) or a refractive index n{z), as depicted in Figure 
|l|. The energy /frequency of the incident particle/wave is 
below the maximum of the potential or cutoff frequency. 
For all experiments we'll consider, the barrier can be 
modelled as a square one, in which V{z) or n(z) is con- 
stant in regions I,II,III but different from one region to 
another. We also assume that V{z) or n{z) is equal in I 
and III and take this value as the reference one. 

The propagation of the particle/ wave through the bar- 
rier is described a by a scalar field ip representing the 



Schrodinger wave function in the particle case or some 
scalar component of the electric or magnetic field in the 
wave case. (The precise meaning of ^ in the case of wave 
propagation depends on the particular phenomenon we 
consider. However, the aim of this paper is to show that a 
common background for all tunnelling phenomena exist). 
Given the formal analogy between the Schrodinger equa- 
tion and the Helmholtz equation ||l2|, this function 
takes the following values in regions I,II,III, respectively: 



= e'*^'" + Re- 



(2) 
(3) 
(4) 



where k and k2 = ix £^re the wave-vectors (p — fik is 
the momentum) in regions I (or III) and II, respectively. 
Note that we have suppressed the time dependent factor 
e*"*. Obviously, the physical field is represented by a 
wave-packet with a given spectrum in lo: 



iPiz,t) = / dtj?7(cj)e*('^^-'^*) 



(5) 



where rj^Lu) is the envelope function. Keeping this in 
mind we use, however, for the sake of simplicity, the sim- 
ple expressions in Eqs. (||), (||), (Q). Furthermore, for the 
moment, we disregard the explicit expression fro k and 
X in terms of the angular frequency uj (or the relation 
between p and E = fiuj) . As well known, the coefficients 
i?, r. A, B can be calculated from the matching condi- 
tions at interfaces: 



(6) 
(7) 



V'/(0) = i/'//(0) , V'//(a) = ipiii{a) 

= ^;,(0) , V;/(a) - ^iiiia) , 

where the prime denotes differentiation with respect to 
z. Substituting Eqs. (|), (|), (|) into (|), (|7|) we are then 
able to find i?, T, A, B and thus the explicit expression for 
the function ip. Here we focus only on the transmission 
coefficient T; its expression is as follows: 



T = [1 - 



(1 - r^) e- 



with: 



X + 
X — ik 



(8) 



(9) 



The interesting limit is that of opaque barriers, in which 
Xa 3> 1. All photonic tunnelling experiments have 
mainly dealt with this case, in which "superluminal" 
propagation is predicted |T6| . Taking this limit into Eq. 
(P) we have: 



1 - i- 



X 



2kx 



(10) 



The quantity (j> in Eq. (^, 
time, is just the phase of T: 



relevant for the tunnelling 
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arctan 



2kx 



(11) 



The explicit evaluation of r in Eq. depends, clearly, 
from the dispersion relations k = k(uj) and x = xi'-^)- 
However, by substituting Eq. (|l^) into (|l|) we are able 
to write: 



= 2 



d k 
duj X 



(12) 



showing that r depends only on the ratio k/x- We can 
also obtain a particularly expressive relation by introduc- 
ing the quantities: 



fci 

Vl 



dk 
doj 



V2 



X 



dx 
duj 



In fact, in this case we get: 



2 

Xk 



X 



ki 



(13) 



(14) 



Note that while ki and fc2 are the real or imaginary wave- 
vectors in regions I (or III) and II, vi and V2 represent 
the "real" or "imaginary" group velocities in the same 
regions. Obviously, an imaginary group velocity (which 
is the case for V2) has no physical meaning, but we stress 
that in the physical expression for the time r in ( [l^ ) only 
the ratio k2/v2 enters, which is a well-defined real quan- 
tity. 

Equations ( p^ ) and ( p^ ) are very general ones (holding in 
the limit of opaque barriers) : they apply to all tunnelling 
phenomena. It is nevertheless clear that peculiarities of a 
given experiment enter in r only through the dispersion 
relations k — k{uj) and x — x('-^) or, better, k(uj)/x(y)- 
As an example of application of the obtained general 
formula, we here consider the case of tunnelling of non 
relativistic electrons with mass m through a potential 
square barrier of height Vq. (In the next sections we then 
study in detail the three types of experiment already per- 
formed). The electron energy \s E = Tiuj (with E <Vo) 
while the momenta involved in the problem are p — hk 
and iq = hk2 — i^X- Iii this case, the dispersion relations 
read as follows: 



k = 



'2mu! 



'2m(Vb - tiLo) 



and thus: 



k 



X y Vq - hw 

By substituting into Eq. ( [T^ ) we immediately find: 
h I 2m 



(15) 
(16) 

(17) 



^/E{Vo - E) n xk 




FIG. 2. Frustrated total internal reflection in a double 
prism. 

III. TOTAL INTERNAL REFLECTION 

The first photonic tunnelling phenomenon we consider is 
that of frustrated total internal reflection [jl^. This is 
a two-dimensional process, but tunnelling proceeds only 
in one direction. With reference to Figure ^, a light 
beam impinges from a dielectric medium (typically a 
prism) with index ni onto a slab with index 712 < ni. 
If the incident angle is greater than the critical value 
9c = arcsinri2/ni, most of the beam is reflected while 
part of it tunnels through the slab and emerges in the 
second dielectric medium with index ni . Note that wave- 
packets propagate along the x direction, while tunnelling 
occurs in the z direction. 

The wave- vectors fci, ^2 in regions I (or III) and II satisfy: 



kl = kl-x^ 



(19) 
(20) 



where kx is the x component of fci or k2 and fc, x £^re as 
defined in the previous section. The dispersion relations 
in regions I (or III) and II are, respectively: 



Kl = — 111 
C 
UJ 

k2 = —n2 

c 



(21) 
(22) 



These equations also define the introduced quantities: 

(23) 
(24) 



c 

Vi = — 

ni 
c 

V2 — 

«2 



It is now very simple to obtain the tunnelling time in the 
opaque barrier limit for this process; in fact, by substi- 
tuting Eqs. (|l|)-(||) into Eq. (|l|) we find: 



Lo xk 



(25) 



(18) Furthermore, using the obvious relations: 



3 



X^^kl sm'0 - fc| 



ril sin^ i 



'■2 I 



we finally get: 



ni sin 



(28) 



(29) 



TT COS 61 y n\ sin ^^ — rij 



This formula can be directly checked with experiments. 
However, we firstly observe the interesting feature of this 
expression which does satisfy the property pointed out 
by Haibel and Nimtz [^. In fact, the time r in Eq. ( |29| ) 
is just given, apart from a numerical factor depending on 
the geometry and construction of the considered experi- 
ment, by the reciprocal of the frequency of the radiation 
used. In a certain sense, the numerical factor can be 
regarded as a "correction" factor to the "universality" 
property of Haibel and Nimtz. 

Several experiments measuring the tunnelling time in the 
considered process have been performed jsj. 
In the experiment carried out by Balcou and Dutriaux , 
two fused silica prisms with ni = 1.403 and an air gap 
(n2 = 1) are used. They employed a gaussian laser beam 
of wave-length 3.39 fim with an incident angle 6 — 45.5°. 
Using these values into Eq. ( |29| ) we predict a tunnelling 
time of 36.8 /s, to be compared with the experimental re- 
sult of about 40 fs. As we can see, the agreement is good 
and the "correction" factor in ( p9| ) is quite important for 
this to occur (compare with the Haibel and Nimtz pre- 
diction of 11.3 fs). 

In the measurements by Mugnai, Ranfagni and Ronchi 
1^, the microwave region is explored, with a signal whose 
frequency is in the range 9 10 GHz. They used two 
parafhn prisms (rii = 1.49) with an air gap (n2 = 1), 
while the incidence angle is about 60°. For this exper- 
iment we predict a tunnelling time of 87.2 ps, while the 
experimental result is 87±7ps 0. 

Finally, we consider the recent experiment performed 
by Haibel and Nimtz Q with a microwave radiation at 
V = S.ibGHz and two perspex prisms (ni = 1.605) sep- 
arated by an air gap (n2 = 1). For an incident angle of 
45°, from ( |2^ ) we predict r = 80.8 ps. The observed ex- 
perimental result is, instead, 117±10ps. In this case, the 
agreement is not very good (while, dropping the "correc- 
tion" factor, Haibel and Nimtz find a better agreement); 
probably this is due to the fact that the condition of 
opaque barrier is not completely fulfilled. 



*Note that the value of 134 ps used by Haibel and Nimtz 
refers to the gap filled with paraffin. In this case no tunnelling 
effect is present. We observe that also for this experiment 
the "correction" factor in ( p9| ) plays a crucial role for the 
tunnelling times 



FIG. 3. A waveguide with an undersized region. 



IV. UNDERSIZED WAVEGUIDE 

Let us now consider propagation through undersized rect- 
angular waveguides as observed in [Q. Also in this 
case, evanescent propagation proceeds along one direc- 
tion (say z) and the results obtained in Sect. II may 
apply. With reference to Figure ^, a signal propagating 
inside a "large" waveguide at a certain point undergoes 
through a "smaller" waveguide for a given distance a. As 
well known [p^ , the signal propagation inside a waveg- 
uide is allowed only for frequencies higher than a typical 
value (cutoff frequency) depending on the geometry of 
the waveguide. In the considered setup, the two differ- 
ently sized waveguides I (or HI) and II have, then, differ- 
ent cutoff frequencies (the first one, loi, is smaller than 
the second one, UJ2), and we consider the propagation of a 
signal whose frequency (or range of frequencies) is larger 
than uji but smaller than uj2. loi < oj < 102- In such a 
case, in the region Q < z < a only evanescent propagation 
is allowed and, thus, the undersized waveguide acts as a 
barrier for the photonic signal. With the same notation 
of Sect. II, the dispersion relations in the large and small 
waveguide are, respectively: 



;fc = 



SO that: 



k 
X 



(30) 
(31) 

(32) 



By substituting this expression into Eq. (|12|), we imme- 
diately find the tunnelling time in the regime of opaque 
barrier (^a 3> 1): 




(33) 



On the contrary to what happens for tunnelling in to- 
tal internal reflection setups, the coefficient of the term 
1/v isn't constant but depends itself on frequency. Thus, 
in the case of undersized waveguides, the assumed "uni- 
versality" property of Haibel and Nimtz cannot apply in 
general; depending on the cutoff frequencies, it is only 
a partial approximate property for frequencies far way 
from the cutoff values (i.e. when the term in the square 
root does not strongly depend on v). 
Let us now compare the prediction ( p3| ) with the ex- 
perimental results obtained in |l|. In the performed 
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FIG. 4. An ideal photonic bandgap device. 



experiment we have microwave radiation along waveg- 
uides whose cutoff frequencies are v\ = 6.56 GHz and 
vi = 9A9GHz, respectively. The radiation frequencies 
are around — 8.7 GHz, so that tunnelling phenom- 
ena occur in the undersized waveguide. By substituting 
these values into Eq. (|3^), we predict a tunnelling time 
of 128 ps, confronting the observed time of about 130 ps. 
As it is evident, also for an undersized waveguide setup 
the theory matches quite well with experiments. Note 
that, despite of the rich frequency dependence in Eq. 
(^), the Haibel and Nimtz property also works quite 
well (although some correction needs), since the central 
frequency value of the radiation used in the experiment 
is far enough from the cutoff values. 



V. PHOTONIC BANDGAP 

The last phenomenon we consider is that of light prop- 
agation through photonic bandgap materials. The ideal 
setup is depicted in Figure ^. Light impinges on a succes- 
sion of thin plane-parallel films composed of N two-layer 
unit cells of thicknesses c?i,d2 and constant, real refrac- 
tive indices ni , n2 , embedded into a medium of index Uq. 
It is known [|l9[ that such a multilayer dielectric mirror 
possesses a (one-dimensional) "photonic bandgap" , that 
is a range of frequencies corresponding to pure imaginary 
values of the wave- vector. In practice, it is the optical 
analog of crystalline solids possessing bandgaps. Increas- 
ing the number of periods will result in an exponential 
increase of the reflectivity, and thus the opaque barrier 
condition can be fulfilled. In general, the study of elec- 
tromagnetic properties of such materials is very compli- 
cated, and the dispersion relation we need to evaluate the 
phase time in the proposed formalism is quite involved 
for physical situations. This study was performed ana- 
lytically in where the dispersion relation (and other 
useful quantities) was derived starting from the complex 
transmission coefficient of the considered barrier. It is, 
then, quite a meaningless issue to get the tunnelling time 
from the dispersion relation obtained from the transmis- 
sion coefficient, while it is easier to have directly the 
phase time r from Eq. (|l|), where </> is the phase of the 
complex transmission coefficient. 



A. Quarter-wave stack 

We first consider the relevant case in which each layer 
is designed so that the optical path is exactly 1/4 of 
some reference wave-length Aq: nidi = ^2^2 = Ao/4. In 
such a case, Aq corresponds to the midgap frequency ujq 
(Aq = 2ttc/luq). This condition is fulfilled in the consid- 
ered experiments Finally, we further assume normal 
incidence of the light on the photonic bandgap material. 
From jl5j we then obtain the following expression for the 
transmission coefficient: 



T [{AC - B) + iADY' , 
where A, B,C, D are real quantities given by: 
sin Nj3 



A = 

B 



sin/3 
sin(iV - 1)(3 
sin/3 



G — a cos — 

. TTUJ 



D 



c sm ■ 



l-^g2 
^02^21^12 

^02^21*12 

2ro2ri2 - - 1 

^02^21^12 



rii + rij 
71,: + rij 



(34) 

(35) 

(36) 
(37) 
(38) 

(39) 
(40) 
(41) 

(42) 
(43) 



sin/3 



^12^21 



'2r?2 



TTUJ 

cos — — 1 



shr — (44) 

0^0 



= 1,2). The phase </> of the transmission coefficient 
thus satisfies: 



tan (h = 



AD 



B - AG 



(45) 



By substituting into Eq. (^, we finally get an analytic 
expression for the tunnelling time of light with frequency 
1/ close to the midgap one i/q for N layers: 



1 1 



sinh N9 



vo 2 sinh(A^ - 1)9 + (6 - a) sinhiVfl 
where is simply obtained from: 



(46) 



5 



sinh^? 



n2 
ni 



"-2 



(47) 



Note that, although the tunnelling behaviour is quite dif- 
ferent if the number of periods N is an even or odd num- 
ber (see, for example, jEffl), the expression for the tun- 
nelling time given in (|46|) (and also in (|4^)) is the same 
in both cases. 

For future reference, we also report the appropriate 
formula for = fc -|- (1/2) (integer k) multilayer di- 
electric mirrors. In practice, this models the case 
of a stratified medium whose structure has the form 
nin2nin2 ■ . ■ nin2ni (note, however, this is an approx- 
imation since, in general, d/2 is not equal to a). In such 
a case, Eq. (E^) is just replaced by: 



1 1 



c cosh NO 



vq 2 cosh(iV -1)0 + {h- a) cosh NO 



(48) 



Let us observe that, similarly to total internal reflection, 
at midgap the time r in Eq. ( p6| ) or (^ is again given 
by the reciprocal of the frequency times a "correction" 
constant factor. 

We now analyze experimental results ||] in the light of 
our theoretical speculations. 

In the experiment performed by Steinberg, Kwiat and 
Chiao, the authors used a quarter-wave multilayer di- 
electric mirror with a {HL)^H structure with a total 
thickness of d = 1.1/im attached on one side of a sub- 
strate and immersed in air. Here, H represents a ti- 
tanium oxide film with rii — 2.22, while L is a fused 
silica layer with n2 — 1.41. Thus, we have approxi- 
mately TV = 5-|-(l/2). As incident light, they employed a 
wave-packet centred at a wave-length Aq = 702nm, corre- 
sponding to the midgap frequency of about A27THz. 
By substituting these numbers in our formula ( |4^ ) we 
predict a tunnelling time r ~ 2.66 /s, corresponding to 
a delay time At, with respect to non tunnelling pho- 
tons propagating at the speed of light the distance d, 
of — 1.01/s. This has to be compared with the exper- 
imental result of At = — (1.47±0.2) /s. However, we 
point out that our analytical prediction is affected by 
two major approximations. The first one is, as already 
remarked, that the experimental sample is not really a 
5 -|- (1/2) periodic structure. A better approximation is 
achieved by using Eq. (|4^) with iV = 6 and subtracting 
the time required for travelling at the speed of light the 
quarter- wave thickness ^2 — Aq /4n2. In this case we have 
T = 2.02 /s or a delay time At = — 1.65/s, which is in 
better agreement with the experimental result. Further- 
more, in our analysis (leading to Eq, (^ ) or (^8|)) there 
is no room for considering an asymmetric structure (like 
the substrate-air one) in which the photonic bandgap ma- 
terial is embedded. This cannot be taken into account in 
an analytic framework, but has to be studied using nu- 
merical matrix transfer method which would give quite a 
good agreement with observations [ p^ . 
Finally, we consider the experiment carried out by Spiel- 
mann et al. ||] on alternated quarter- wave layers of fused 



silica L and titanium dioxide H having the structure 
of (substrate) (HL)'' (air) with N = 3,5,7,9,11. They 
used optical pulses of frequency 375THz corresponding 
to the midgap frequency of their photonic bandgap ma- 
terial. Obviously, increasing N wc have a better re- 
alization of opaque barrier condition. From Eq. (|4^ ) 
with TV = 11 (note, however, that for > 5 the fac- 
tor sinh(Af — 1)0/ sinh A^6' is almost constant) we have a 
tunnelling time of 2.98 /s to be compared with the ob- 
served value of about 2.71 fs. We address the fact that, 
apart the presence of the asymmetric substrate-air struc- 
ture which introduces some approximation as discussed 
above, in the considered experiment the incidence of the 
light on the sample is not normal, being ?» 20° the angle 
between the axis of the sample and the beam propaga- 
tion direction. In this case, the described computations 
are only approximated ones and, again, the exact result 
can be obtained only through numerical implementation. 
Nevertheless, also within the limits of our calculations, 
the agreement between theory and experiment is quite 
good. 

A final comment regards the predictions of the "univer- 
sality" property proposed by Haibel and Nimtz. Neglect- 
ing the "correction" factor in Eq. ( ^ ) would yield the 
values of At = — 1.33/s and t = 2.67 /s for the de- 
lay time in the Steinberg, Kwiat and Chiao experiment 
and the transit time for the Spielmann et al. experi- 
ment, respectively. In both cases, the agreement with 
the observed values seems better than our approximated 
predictions, showing that the presence of an asymmetric 
substrate-air structure (and the non normal incidence in 
the second experiment) pushes up the "correction" factor 
inEq. (El). 



B. Distributed Bragg Reflector 

We now relax the assumption of a quarter-wave stack 
nidi = ^.2^2 = Ao/4 but, for simplicity, we consider only 
the case in which the photonic bandgap structure is em- 
bedded into a material whose refractive index hq is equal 
to that of one of the two layers of the periodic structure, 
that is no — n2. We again assume normal incidence of 
the light on the photonic bandgap material. In this case 
the transmission coefficient T and its phase (j) have the 
expressions as in ( p^ ) and (^5|), where A,B are given by 
(H), m and 



C = a cos TT LU — bcOSTTQ^UJ 

D — — a siiiTT Lu + b siu-TT Q-Uj 



^4 



sin/3 = 



^12^2] 



nidi ± n2d2 



+ Q + R 



(49) 
(50) 

(51) 
(52) 



6 



P = sin TT f7_ bj 



4 „;„2 
2 

Q = 2r^2 (cos TT r2-|_ w cos tt r2_ u 
R = sin^ TT ri4- w 



(53) 
(54) 
(55) 



By substituting into Eq. (^ we obtain the tunnelling 
time relative to an A'^-layer structure: 



I X -Y 

^ = 

X^F sin^ (5 cos iV/3 sin N (i 

Y = G (cos /3 cos NI3 sin N^-N sin /3) 

Z = 2 sin/3 (1)2 sin^ A^/? + sin^ /3 cos^ iV/3) 

F = a r2_|_ cos 7ri7_|_cj — b lu cos 7ri7_ti; 

G ^ a? fi+ cj sin^ vrfi i + b^ w sin^ 7rr2_L 



(56) 

(57) 
(58) 
(59) 
(60) 



— 2a6 (il+ + il_) ti) sinTril+ii; sin^ 7rf2_CLi (61) 

Note that, again, the formula above for r holds both for 
even N and for odd N . 

The obtained expression for the tunnelling time can be 
directly tested by analyzing the experiment carried out 
by Mojahedi, Schamiloglu, Hegeler and Malloy [^. In 
this experiment the authors used a (ID) photonic crys- 
tal composed of 5 polycarbonate sheets with refractive 
index ni = 1.66 and thickness di = 1.27 cm separated by 
regions of air n2 = I with thickness d2 — 4.1cm. The 
bandgap was tuned to the main frequency component 
{i^ = 9.68 GHz) of the incident microwave pulse. By 
measuring both the signal travelling through the pho- 
tonic bandgap structure and the one propagating in free 
space, the authors found that the pulse undergoing tim- 
nelling has a delay time At = — (440±20)ps with respect 
to the other signal. By using Eq. ( ^6|) with the above 
numbers we predict a tunnelling time of 320 [], corre- 
sponding to a delay time of At = — 438ps, which is in 
excellent agreement with the reported experimental re- 
sult. 

We point out that, in this case, the simple l/v law pro- 
posed by Haibel and Nimtz does not work, since it would 
predict a tunnelling time r = 103ps or At = —655ps. 
This can be easily explained by looking at Eq. (p6|). In 
fact, we immediately recognize that the "correction" fac- 
tor in this equation is strongly frequency-dependent and, 
for the frequency of the light used in the considered ex- 
periment, it is sensibly bigger than one. 



VI. CONCLUSIONS 

In this paper we have scrutinized the recently proposed 
H "universality" property of the photonic tunnelling 



^Such a result was also obtained in [g] using a formalism 
described in |4[ which is different from the one proposed here. 



time, according to which the barrier traversal time for 
photons propagating through an evanescent region is ap- 
proximately given by the reciprocal of the photonic fre- 
quency, irrespective of the particular setup employed. To 
this end, the transit time in the relevant region, defined 
here as in Eq. (Q), needs to be computed for the different 
explored phenomena, and in Sect. II we have given gen- 
eral expressions for this time in the opaque barrier limit. 
The peculiarities of a given photonic setup enter in these 
expression only through the dispersion relation relating 
the wave-vector and the frequency. More in detail, we 
have shown how the knowledge of the ratio between the 
wave- vectors in the barrier region and outside it, as a 
function of the photon frequency, is sufficient to evaluate 
the transit time r in Eq. (|l2|). 

Several specific cases, corresponding to the different 
classes of experimentally investigated phenomena, have 
then been considered. In particular, in Sect. Ill we have 
studied light propagation in a setup in which the evanes- 
cent region is provided by total internal reflection, while 
in Sect. IV the propagation through undersized waveg- 
uides has been considered and, finally, in Sect. V the case 
of a photonic bandgap has been analyzed. The relevant 
results for the three mentioned phenomena are given in 
Eqs. (H), (HI) and (|^) or respectively. As can be 
easily seen from these expressions, the frequency depen- 
dence of the tunnelling time for the cases of total internal 
reflection and quarter-wave photonic bandgap is just as 
predicted by the property outlined by Haibel and Nimtz 
Ig, although we have derived a "correction" factor de- 
pending on the geometry and on the intrinsic properties 
of the sample (this factor is not far from unity). On 
the contrary, such a factor is frequency dependent for 
undersized waveguides and distributed Bragg reflectors, 
revealing a more rich dependence of r on than the sim- 
ple l/v one (see Eq. (p3|)). We can then conclude that 
the "universality" property of Haibel and Nimtz is only 
an approximation, but it gives the right order of magni- 
tude for the tunnelling time. This conclusion holds also 
for undersized waveguide propagation, provided that the 
photon frequency is far enough from the cutoff frequen- 
cies. 

We have then calculated the tunnelling times for the dif- 
ferent existing experiments and compared the theoretical 
values with the observed ones. Results are summarized in 
Table 1, where we also report the Haibel and Nimtz pre- 
diction 1/i^. From these we can see that, in general, the 
agreement of our prediction with the experimental values 
is satisfactory. As pointed out in the previous section, the 
calculations performed here for photonic bandgap mate- 
rials assume some approximations in treating the com- 
plex sample, which are nevertheless required to obtain 
analytical expressions. Our prediction then suffer of this 
and, in the case in which the setup is designed to ver- 
ify the quarter-wave condition nidi = 712^2 = Ao/4, the 
simple l/v rule fits better with experiments while, for 
general photonic bandgap structures, the tunnelling time 
displays a very complicated dependence on frequency. In 
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TABLE I. Comparison between predicted and observed 
tunnelling times for several experiment (FTIR, UWG and 
PBG stands for frustrated total internal reflection, undersized 
waveguide and photonic bandgap, respectively). Te^p is the 
experimental result while Tth is our prediction as from Eqs. 
(p9|), (js^) and (^) or (p^. For reference to the Haibel and 
Nimtz property, we also report the value 



Phenomenon Experiment 






Tth 




FTIR 


Balcou et al. [ 




11.3 /s 


36.8 fs 


~ 40 /s 


FTIR 


Mugnai et al J 




100 ps 


87.2 ps 


87±7ps 


FTIR 


Haibel et al. f 




120 ps 


81 ps 


117±10ps 


UWG 


Enders et al. | 




115ps 


128 ps 


~ 130 ps 


PBG(Ao/4) 


Steinberg et al. 


L 2.34 /s 


2.02 fs 


2.20±0.2/s 


PBG(Ao/4) 


Spielmann et al. 


1 2.67 /s 


2.98 /s 


~ 2.71 fs 


PBG 


Mojahedi et al. [F 


1 103 ps 


320 ps 


318±20 ps 



this last case, as well as in all other non photonic bandgap 
experiments, the "correction" factor introduced in this 
paper is quite relevant for the agreement with observa- 
tions to be good. 
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